The Spintronic Properties of Rare Earth Nitrides 
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The electronic structure of the rare earth nitrides is studied systematically using the ab-initio 
self-interaction corrected local-spin-density approximation (SIC-LSD). This approach allows both 
a localised description of the rare earth /—electrons and an itinerant description of the valence 
electrons. Localising different numbers of /-electrons on the rare earth atom corresponds to different 
valencies, and the total energies can be compared, providing a first-principles description of valence. 
CeN is found to be tetravalent while the remaining rare earth nitrides are found to be trivalent. We 
show that these materials have a broad range of electronic properties including forming a new class 
of half-metallic magnets with high magnetic moments and are strong candidates for applications in 
spintronic and spin-filtering devices. 
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I. INTRODUCTION 

Rare earth (RE) materials have high magnetic mo- 
ments and form a wide range of magnetic structures. It 
is the occupancy of the highly localised 4/-states that 
determines the magnetic properties while the other elec- 
tronic properties are determined principally by the itiner- 
ant s — d electronsi The standard method of determining 
the electronic properties of magnetic materials is density 
functional theory with a local spin density approximation 
(LSD) for the exchange-correlation energy.^ Tradition- 
ally this method has been applied to the s — d electrons 
in rare earths while the /-electrons have been treated in 
an atomic model. Attempts to include the 4/-electrons 
within the LSD have had only very limited success be- 
cause the LSD is insufRcient to describe the strong corre- 
lations experienced by the rare earth /-electrons. In re- 
cent years more advanced methods of electronic structure 
determination such as LSD plus Self-Interaction Correc- 
tions (SIC)^ and LDA-HL^ methods have sought to rem- 
edy this problem and have met with considerable success. 

Currently there is a great deal of interest in spintron- 
ics and spin-filtering from both the fundamental and ap- 
plied point of view. There has been speculation in the 
literature that rare earth nitrides may form half-metallic 
ferromagnetsiSi^ This is surprising because, based on a 
simple ionic model, trivalent rare earth nitrides would 
be expected to be insulators with a similar electronic 
structure to the divalent rare earth chalcogenides. How- 
ever if it is the case that they are half-metallic then they 
have potential applications in spintronic devices. To our 
knowledge, there has been no systematic study of the 
electronic structure of these materials at all. This is par- 
ticularly surprising given that they all crystallise in the 
simple sodium chloride structure.^ 

It is the purpose of this paper to discuss our calcu- 
lations of the electronic structure of rare earth nitrides 
and thus to gain a deeper understanding of these materi- 
als and their potential applications. We have performed 
SIC-LSD calculations for all the rare earth nitrides, in 
the ferromagnetic state, in both the divalent and triva- 



lent states. Thus, in this systematic study of rare earth 
nitrides we shall be concerned with the nominal valence 
of the rare earth ion as we progress through the series. 
Will these materials turn out to be all trivalent? What 
is the character of the rare earth / and d hybridization 
with the nitrogen p states? Obviously, the position of the 
empty /-states will be an important factor in determin- 
ing this. By performing this systematic study, we will be 
able to determine the trends. 

The remainder of the paper is organized as follows. 
After a brief summary of our calculational method in 
section II, in section III we discuss the detailed electronic 
structure of these materials, and trends in behaviour as 
one proceeds across the Periodic Table. We conclude the 
paper in section IV, indicating what our results imply for 
device applications. 



II. SIC-LSD 

The SIC-LSD approximation'^ is an ab-initio electronic 
structure scheme, that is capable of describing localiza- 
tion phenomena in solids AiSiil In this scheme the spu- 
rious self-interaction of each occupied state, which is in- 
herent in any local approximation to Density Functional 
Theory (DFT)ji2ii^ is subtracted from the conventional 
LSD approximation to the total energy functional 
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where a labels the occupied states and S^^^ is the 
self-interaction correction for state a. This leads to a 
greatly improved description of static Coulomb correla- 
tion effects over the LSD approximation. Demonstra- 
tions of the advantages of the SIC approach have been 
made in studies of the Hubbard model j^^i^^ in appli- 
cations to 3d oxides such as monoxides jiSiii high tem- 
perature cuprateSfiS*iii2 and magnetiteji^ to /-electron 
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and to solid hydrogen 
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In the SIC-LSD method it is necessary to minimise the 
total energy with respect to the number of localised elec- 
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trons which also leads to a determination of the nominal 
valence defined as the integer number of electrons avail- 
able for band formation 
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where Z is the atomic number of the rare earth and N^ore 
is the number of atomic core electrons. Nsic is the num- 
ber of /-states for which the self-interaction correction 
has been removed and is determined so that equals 3 
for trivalent, and 2 for divalent, systems. 

Further details of the SIC-LSD method are discussed 
elsewhere!^ 



III. RESULTS AND DISCUSSION 




A. Energetic Properties 

We have performed calculations for all the rare earth 
nitrides in ferromagnetic ordering and both the divalent 
/"+^ and trivalent /" configurations (and for CeN the 
tetravalent /""""^ configuration), and hence have deter- 
mined the ground state energy and valence from first 
principles. We have found that the ground state for Ce 
in CeN is tetravalent while all the other rare earths are 
trivalent in the nitrides. In Fig. ^ we show the differ- 
ence in energy between the divalent and trivalent states 
of all the rare earth nitrides. This follows the expected 
trend of strongly trivalent state at the start of the series 
in both spin channels, and a decreasing energy difference 
between the trivalent and divalent state as the filling of 
the /-shell occurs. However, none of the nitrides attains 
the divalent state. Furthermore the number of occupied 
itinerant /-electron states stays well below the 0.7 re- 
quired for the materials to become divalent, as observed 
by Strange et al.,'^^ because of the strong hybridisation 
with the nitrogen p states. The lattice constants for the 
/ configuration that yields the lowest total energy can be 
evaluated by looking for the minimum in total energy as 
a function of lattice size and these are shown in Fig. |21 
along with the experimental valuesp2& The jump in lattice 
constant between Ce and Pr is due to the valence change 
from tetravalent to trivalent. For Ce, Pr and the heavy 
rare earths, theory and experiment agree very well both 
for the absolute values and the trend. For the remain- 
ing light rare earths the absolute agreement is still very 
good but there is a trend apparent in the theory that is 
not clear in the experimental results. The most proba- 
ble cause of this is experimental uncertainty which can 
arise in these materials very easily, particularly as many 
of these measurements were performed at a time when 
samples of rare earth materials were known to have very 
limited purity»2Si 
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FIG. 1: The difference in total energy between the divalent 
and trivalent states of the rare earth nitrides. The differences 
are all positive indicating that the trivalent state is preferred 
over the divalent state in all cases in agreement with the pub- 
lished literature. Unlike in the previous publicatiouE^'' no 
spin-orbit coupling has been included in the present calcula- 
tions. 
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FIG. 2: The lattice constant for the rare earth nitrides, crosses 
are the calculated values and the squares are experimentally 
determined for thin films of the materialp— 
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TABLE I: Various components of the spin magnetic moment 
(Ms) of the rare earth nitrides. All values are in Bohr mag- 
netons (ps). Where the numbers do not quite add up to 
the total the remaining spin moment is in the empty spheres 
used in the calculation. The rare earth orbital moment (Ml) 
assumes that the rare earth ions obey Hund's rules. 



Material 




Ms 




Mi 




RE 


N 


Total 


RE 


CeN 


0.0 


0.0 


0.0 


0.0 


PrN 


2.07 


-0.08 


2.00 


-5.0 


NdN 


3.10 


-0.11 


3.00 


-6.0 


PmN 


4.13 


-0.14 


4.00 


-6.0 


SmN 


5.22 


-0.24 


5.00 


-5.0 


EuN 


6.30 


-0.30 


6.00 


-3.0 


GdN 


7.01 


-0.04 


7.00 


-0.0 


TbN 


5.97 


0.01 


6.00 


3.0 


DyN 


4.93 


0.05 


5.00 


5.0 


HoN 


3.91 


0.08 


4.00 


6.0 


ErN 


2.90 


0.09 


2.99 


6.0 


TmN 


1.83 


0.12 


1.96 


5.0 


YbN 


0.79 


0.14 


0.94 


3.0 



B. Magnetic Properties 

Now let us look at how the magnetism varies as we pro- 
ceed along the rare earth series. The total and species- 
decomposed spin magnetic moments are displayed in Ta- 
ble HJ There we also present the rare earth's orbital mo- 
ments. It should be noted that there is some arbitrari- 
ness involved in the values of the spin magnetic moments 
due to the Atomic Sphere Approximation. However, we 
found this to be small and moreover to make comparisons 
between different systems consistent we kept the ratio of 
the rare earth/nitrogen sphere volumes constant. 

With the exceptions of ErN, TmN and YbN, the spin 
magnetic moments of the rare earth nitrides take on an 
integer value. This indicates that these systems are either 
insulating (Tb-, Dy-, and Ho- nitrides) or half- metallic 
(Pr- to Gd-compounds) . CeN is a non-magnetic metal, 
and the last three compounds of the series are metallic 
in both spin-channels. These results are as one would 
expect; the spin magnetic moment is dominated by the 
rare earth /-electrons, with some hybridisation yielding 
small contributions from the rare earth s—d electrons and 
the nitrogen p states. This indicates that the nitrogen p 
states occur in the same energy range as the valence rare 
earth states, allowing hybridization to occur. It is inter- 
esting to note that the contribution from the nitrogen 
atom changes sign half way through the series, in which 
it follows the RE's orbital moment. It appears that the 
nitrogen moment wants to point antiparallel to the par- 
tially occupied /-spin channel. An explanation of this 
hybridisation phenomenon will be given in sections HI D 
and E. 
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FIG. 3: The density of states for Terbium Nitride in the triva- 
lent state, (Tb Site = full line. Nitrogen site = dashed line). 



C. Density of States 

In Fig. Owe show typical densities of states of RE ni- 
trides to illustrate the main features. Around -1.5 Ryd 
are the rare earth 5p bands. Above these, just above -1.0 
Ryd, are the nitrogen 2s bands. Just below the Fermi 
energy Ef, we reach bands which are predominantly ni- 
trogen p-like, but with a substantial hybridisation with 
rare earth s — d and /-states. There is a small gap or 
mimimum around the Fermi energy. Above Ef are the 
mainly rare earth s — d bands. Superimposed on this are 
the rare earth /-bands. The fully occupied majority spin 
/-bands fall below the nitrogen 2s bands and push the 
rare earth bp bands lower in energy. The occupied mi- 
nority spin /-bands sit just above the nitrogen 2s bands 
and there is essentially a single peak above E / which are 
the unoccupied /-states. Below half filling the occupied 
majority spin peak moves to just above the nitrogen 2s 
bands, and just above Ef there are two /-peaks, firstly 
the unoccupied majority spin bands and then the com- 
pletely unoccupied minority spin /-bands. 

Let us look more closely at the electronic properties 
of the rare earth nitrides. Table Hll shows the electronic 
properties of each material at the Fermi energy. It is 
clear from the table that most of the light rare earth ni- 
trides are found to be half- metallic. Only CeN is metallic 
because it exists in the tetravalent state (in fact the triva- 
lent state is also just metallic). TbN, DyN and HoN are 
found to be narrow gap insulators and ErN, TmN and 
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TABLE II: Spin resolved band gaps and densities of 
states at the Fermi energy for the rare earth nitrides. 
Band gaps are in Rydbergs and densities of states are in 
states/Rydberg/formula unit. 



Material Band Gap DOS DOS 

Spin-up Spin-down Spin-up Spin-down 



CeN 








15.37 


15.37 


PrN 





0.039 


0.0001 
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PmN 





0.076 


31.77 





SmN 





0.095 


154.56 





EuN 





0.107 


69.85 





GdN 





0.082 


0.065 





TbN 


0.008 


0.052 








DyN 


0.018 


0.058 








HoN 


0.031 


0.004 








ErN 








0.682 


69.57 


TmN 








1.52 


220.78 


YbN 








1.72 


93.18 



TABLE III: The spin resolved and ^-decomposed densities 
of states at the Fermi energy for the rare earth nitrides in 
states/Rydberg/formula unit. 



Material 




DOS-up 




DOS-down 






RE s-d 


RE f 


N p 


RE s-d 


RE f 


N p 


CeN 


2.10 


11.19 


0.13 


2.10 


11.18 


0.13 


PrN 


0.0001 

















NdN 


0.020 


0.003 


0.030 











PmN 


0.061 


31.02 


0.40 











SmN 


0.44 


151.42 


2.06 











EuN 


0.64 


64.76 


3.24 











GdN 


0.013 


0.002 


0.040 











TbN 




















DyN 




















HoN 




















ErN 


0.017 


0.037 


0.605 


0.069 


68.45 


0.634 


TmN 


0.029 


0.070 


1.30 


0.28 


218.81 


1.26 


YbN 


0.098 


0.08 


1.46 


0.37 


89.71 


2.10 


YbN are 


metallic 


in both 


spin- 


channels. 







In Table ITTll we show the key spin- and /—decomposed 
densities of states at the Fermi energy. This table shows 
that the principal character of the electrons a,t Ej is rare 
earth s ~ d and nitrogen p-like in PrN, NdN, and GdN, 
while it is dominated by the majority spin rare earth 
/-electrons in PmN, SmN, EuN, and the minority spin 
rare earth /-electrons in ErN, TmN and YbN. When the 
/ contribution is large we also note a relatively large N 
p component, reflecting the strong N p - rare earth / 
hybridization. From the results above it is also clear that 



there is a wide range of electronic properties in the rare 
earth nitrides which we discuss in detail in the following 
section. 

To discuss the properties of these materials more fully 
it is necessary to understand the details of the density 
of states around Ef. For this reason we show in Fig. 0] 
the spin-resolved densities of states on each site for SmN, 
DyN and TmN. These materials were chosen as typical 
examples of the half-metallic, insulating and metallic rare 
earth nitrides respectively. 



D. The Light Rare Earth Nitrides 

We now explain how the rare earth nitrides can ex- 
hibit such a wide range of electronic properties in a se- 
ries of materials that ostensibly have similar outer elec- 
tronic structures. In previous papers^SiSi we have calcu- 
lated the electronic structure of rare earth chalcogenides. 
In agreement with observation and with the naive ionic 
model, we found that they are all divalent insulators. 
One might have expected the trivalent nitrides to be sim- 
ilarly insulating, but, as Table ITTl shows, this is not the 
case. After the SIC-LSD has been applied to the rare 
earth nitrides the three valence electrons of the trivalent 
rare earth can fill the three holes of the N p band and 
create a full p band. This would lead to insulating be- 
haviour in a similar manner that the two electrons of the 
divalent rare earth in the rare earth chalcogenides fill the 
two holes in the chalcogenide p band. In the nitrides 
there are three nitrogen p electrons, in the chalcogenides 
there are four p electrons. The difference in the number 
of p electrons means that the p bands sit closer to the 
Fermi energy in the nitrides than in the chalcogenides. 
In trivalent CeN the nitrogen p bands are close enough 
to Ef to hybridise slightly with the low energy tail of the 
rare earth s-d bands, causing the materials to be metal- 
lic. As we proceed across the Periodic Table from CeN 
to GdN we fill the majority spin /-states and create an 
exchange field that is felt by the other electrons in the 
material. In the rare earth chalcogenides this has lit- 
tle effect because the filled bands are well down in the 
potential well. However, in the nitrides its effect is more 
pronounced and leads to a significant spin-splitting of the 
nitrogen p-bands. 

The SIC splits the / band manifold into two, local- 
ized and band-like / electrons. The SIC /-bands are 
completely localised and are unable to hybridise signifi- 
cantly. The non-SIC / states have a degree of itineracy 
that allows hybridisation. We find that N p-rare earth 
/ hybridised states play a role in the electronic bonding. 
This is clear in Fig. 01 In a(ii) panel for example, the N 
p-states, which lie within 0.2 Ry of Ef, lead to structure, 
in this energy range, in the Sm density of states. This 
is an energy window where / states do not occur, a p-f 
hybrised bonding state has been formed. This occurs to 
some extent in all the rare earth nitrides. 

Of course, there is also &p-f antibonding state and it is 
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(E - Ef) (Rydbergs) (E - Ef) (Rydbergs) 



b(I) 




(E - Ef) (Rydbergs) 




(E - Ef) (Rydbergs) 



c(i) 





(E - Ef) (Rydbergs) 



(E - Ef) (Rydbergs) 



FIG. 4: The density of states around the Fermi energy for (a) SmN, (b) DyN, and (c) TmN. The (i) refers to the rare earth 
site and the (ii) to the nitrogen site. Full lines represent majority spin and dashed lines are minority spin 



this that is key to understanding the trends in electronic 
properties in this series of materials. The upper tail of 
the majority spin p-bands rises to above the Fermi energy 
and hybridises with the empty majority spin bands just 
above Ef to form the antibonding state. For SmN this 



can be seen clearly in a{i) and a{ii) in Fig. ^ Just above 
Ef the majority spin nitrogen p density of states in a{ii) 
has a small peak that corresponds with the rare earth 
/-bands in a{i). 

The exchange field keeps the minority spin /-bands 
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well above Ef and hybridisation of them with the nitro- 
gen p-bands is very small. In PrN and NdN the empty 
majority spin /-bands are well above Ef and so the elec- 
tronic properties are dominated by the majority spin rare 
earth s-d bands. In PmN, SmN and EuN the empty ma- 
jority spin /-bands have lowered and are very close to 
the Fermi energy. The density of states at Ef is then 
dominated by p-f bands. 

In the minority spin channel there is no significant p — 
f hybridisation and so the predominantly p-bands are 
entirely filled and the s-d bands are empty. There is an 
energy gap at the Fermi energy. The minority spin bands 
behave much more as expected in the naive ionic model. 

A consequence of this behaviour is that the nitrogen 
majority spin bands are not fully occupied while the mi- 
nority spin bands arc. This means that the nitrogen mo- 
ment is antiparallel to the rare earth moment. 



E. The Heavy Rare Earth Nitrides 

In the heavy rare earth nitrides a repeat of the be- 
haviour of the light rare earths might be expected with 
the minority spin bands replacing the majority spin ones. 
This does indeed occur. Inspection of Figs. EJd and 4c 
verifies that the comparable features can be seen in the 
minority spin densities of states. However, other effects 
also come into play in the heavy rare earth nitrides that 
modify the electronic properties significantly. 

In GdN, the majority spin /-bands are all filled and 
the minority spin bands are well above E f . The exchange 
field continues to split the nitrogen p-band. The majority 
spin p-band hybridises with the Gd s — d bands. There 
is a minimum in the density of states very close to Ef, 
but it does not quite reach zero. This is very similar to 
what happens in CeN 

In TbN, the minority spin /-bands begin to fill up. 
This pushes the majority spin /-bands closer to the nu- 
cleus and so increases the efficiency of the screening of 
the outer states so they are less strongly bound to the 
nucleus. Because of this the rare earth s — d states are 
higher above the Fermi energy in the heavy rare earths 
than in the light rare earths. In TbN, DyN and HoN 
then, there are no bands at the Fermi energy and there- 
fore they are small gap insulators. In ErN, TmN and 
YbN the minority spin /-peak has come down to close to 
Ef and there is a bonding and an antibonding minority 
spin state as in the light rare earths. As the /-peak comes 
down the antibonding state gets very close to Ef. The 
Fermi energy enters the major minority spin /-peak. To 
compensate for this some hybridised majority spin p ~ d 
character has to rise above Ef creating a non-zero density 
of states in both spin channels. 

The p — f hybridisation pulls some minority spin char- 
acter above the Fermi energy and so for the heavy rare 
earths the nitrogen moment is parallel to the rare earth 
moment. 

For most rare earth nitrides very little experimental 



data appears to be available. However, for YbN several 
papers have appeared in the literature reporting experi- 
ments which can be compared with our work. Trivalency 
has been established from the magnetic moment i2i Ott 
et al., Sakon et al. and Takeda et al?^ reported mea- 
surements consistent with heavy fermion behaviour in 
YbN. However it was shown later by Monnier et alM- 
that these results were also consistent with a scenario 
in which there is interplay between the crystal field and 
the Kondo effect for an isolated magnetic impurity. The 
SIC-LDA method does not include the dynamical effects 
necessary to describe the Kondo effect or heavy fermion 
behaviour, so we are unable to comment on this point 
directly. However, we can point out that our results do 
exhibit the p-f mixing necessary for the p-f Kondo ef- 
fect that explains the low temperature anomalies in the 
magnetic susceptibility and specific heat^ 

The low temperature magnetic structure of YbN was 
established to be antiferromagnetic HI by Donni et a/.— 
Spectroscopic investigations of YbN include core level x- 
ray photoemission spectroscopy (XPS)"^^ which revealed 
no evidence of valence mixing, and optical spectroscopy^^ 
which was able to locate the empty / state 0.2 eV above 
the Fermi energy. They also found the bare occupied /- 
states to be around 6.5 eV below the Fermi energy, which 
is not in good agreement with our results, although there 
are always considerable uncertainties in locating ground 
state bands using spectroscopic techniquesi^ The XPS 
experiments also showed that samples exhibiting a small 
amount of non-stoichiometry can show strongly enhanced 
Kondo features. The experimental situation for YbN has 
been well summed up by Wachter— who shows that small 
amounts of doping could raise the Fermi energy to the 
empty /-state in YbN and yield intermediate valence or 
heavy fermion behaviour. In our calculations the posi- 
tion of the empty /-state is in close agreement with the 
optical spectroscopy results and so we are also able to 
suggest that. Furthermore the empty /-states in TmN 
are actually at the Fermi energy as indicated in Table HI. 
Extrapolating from the results for YbN suggests that un- 
doped TmN might exhibit heavy fermion or intermediate 
valent behaviour as well. 



IV. CONCLUSIONS 

We have seen that the rare earth nitrides display a 
wide range of electronic properties despite having the 
same crystal structure with only a small variation in lat- 
tice constant and superficially similar electronic struc- 
tures. They show insulating (semiconducting), half- 
metallic and full metallic behaviours. This is a conse- 
quence of the N p and the band-like rare earth / states 
occuring in the same energy window, in the vicinity of the 
Fermi level. This leads, at E/, to strong hybridization of 
these states. With the high magnetic moments and the 
high densities of states in one spin channel that these ma- 
terials exhibit, it suggests that alloying may enable us to 
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fabricate materials with a wide and continuous range of and spin filtering applications, 
useful properties, particularly with regard to spintronic 
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